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ABSTRACT 



\ In this work the double vector meson production in ultraperipheral heavy ion collisions is 

investigated within the BFKL approach. The integrated cross sections and event rates for the 
O ■ processes AA — > V1V2 AA (Vi = p, lj, 4>, J/^, T) are computed and theoretical estimates for 
scattering on both light and heavy nuclei are given for energies of RHIC and LHC. 
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43 • 1 Introduction 

> : 

Photon collisions offer unique possibility to probe QCD in its high energy limit [1]. The 
simplicity of the initial state and the possibility of studing several different combinations of 
final states make this process very useful for studying the QCD dynamics in the limit of high 
center-of-mass energy y^i and fixed momentum transfer t. In principle, in this regime the 
QCD dynamics is determined by the QCD Pomeron, with the evolution described by the 
BFKL equation [2]. One typical process where the QCD Pomeron can be tested is the vector 
meson pair production in 77 collisions [3]. For heavy vector mesons, this cross section can 
be calculated using the perturbative QCD methods. First calculations considering the Born 
two-gluon approximation have been done in Refs. [3] and after considering the solution of 
the BFKL equation [4]. Recently, we have studied the double vector meson production in 77 
collisions, where the the color singlet t-channel exchange carries large momentum transfer [5]. 
The presence of a large momentum transfer allow us to calculate the light and heavy vector 
mesons production cross sections for real photon interactions. We have considered the non- 
forward solution of the BFKL equation at high energy and large momentum transfer [6] and 
estimated the total and differential cross section for the process 77 — > V1V2, where V\ and V2 
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can be any two vector mesons (V^ = p,u,<f), J/\I/,Y). Our results indicate that this process 
could be used in the future 77 colliders to constrain the QCD dynamics. 



On the other hand, in Refs. [7, 8, 9, 10] we have shown that several processes which can 
be observed in ultraperipheral heavy ion collisions (UPC's) at RHIC and LHC can be used to 
constrain the QCD dynamics. In these collisions the ions do not interact directly with each 
other and move essentially undisturbed along the beam direction. The only possible interaction 
is due to the long range electromagnetic interaction and diffractive processes (For a review see, 
e. g. Refs. [11,12, 13]). Consequently, it is possible to study photon- nucleus and/or two-photon 
interactions in these coherent collisions. An important aspect which should be emphasized is 
that relativistic heavy-ion collisions are a potentially prolific source of 77 collisions at high 
energy colliders. The advantage of using heavy ions is that the cross sections varies as Z 4 a 4 
rather just as a 4 . Moreover, the maximum 77 collision energy W 77 is about 6 GeV at 

RHIC and 200 GeV at LHC, where Ra is the nuclear radius and 7 is the center-of-mass system 
Lorentz factor of each ion. In particular, the LHC will have a significant energy and luminosity 
reach beyond LEP2, and could be a bridge to 77 collisions at a future e + e~ linear collider. For 
two-photon collisions, the cross section for the reaction AA — > AA V 1 V 2 will be given by (See 



where the photon energy distribution n{uj) is calculated within the Weizsacker- Williams ap- 
proximation [11]. In general, the total cross section AA — > AA77 — > AAX, where X is the 
system produced within the rapidity gap, factorizes into the photon-photon luminosity dC^/dr 
and the cross section of the 77 interaction, 



where r = s/s, s = W 2 is the square of the center of mass (c.m.s.) system energy of the 
two photons and s of the ion-ion system. The 77 luminosity, d £ 77 (r) jdr , is given by the 
convolution of the photon fluxes from two ultrarelativistic nuclei. Here, we consider the photon 
distribution of Ref. [14], providing a photon distribution which is not factorizable. The authors 
of [14] produced practical parametric expressions for the differential luminosity by adjusting the 
theoretical results. The comparison with the complete form is consistent within a few percents. 
The approach given above excludes possible final state interactions of the produced particles 
with the colliding particles, allowing reliable calculations of UPC's. 

In order to estimate the double meson production in UPC's it is only necessary to consider 
a suitable QCD model for the subprocess at photon level. In Refs. [7, 15] we have analyzed 
several distinct scenarios for the double J/^, pJ/^ and double p production, respectively. In 
particular, in Ref. [7] we have estimated the double J/^/ production considering the solution of 
the BFKL equation at zero momentum transfer and using a small t approximation in order to 
calculate the total cross section at photon level. Moreover, we have evaluated the impact of non- 
leading corrections to the BFKL equation and compared our predictions with those obtained 
at Born two-gluon approximation. On the other hand, in Ref. [15] we have estimated the 



Fig. 1) 



VAA-^AAVx v 2 — 



/ — M^i) / — ^2(^2) 



(1) 




(2) 



2 



Figure 1: The process AA — > V1V2AA and the corresponding kinematical variables. The two- 
photon scattering is sketched, showing the photon-meson impact factors T lVl < 2 and the BFKL 
Green's function /Cbfkl- The grey blobs represent the colliding nuclei. 

pJ/ty production considering the approach proposed in Ref. [16], where the 77 — > pJ/^ cross 
section is given in terms of the gluon distribution on the light meson. Moreover, the double pp 
production has been estimated using the pomeron-exchange factorization theorem [17], which 
allows relate the 77 — > pp cross section with those for J/^/J/^/ and pJ/^ production. Our 
main result in those analyzes was to show that the QCD dynamics implies a large enhancement 
of the cross sections, which become these processes observable in future AA colliders and the 
ultraperipheral heavy ion collisions an important alternative to study the QCD dynamics at 
high energies. It strongly motivates a more detailed study of the double meson production in 
UPC's considering more precise theoretical formalisms for the 77 — > V1V2 process, as well as to 
obtain predictions for other combinations of vector mesons in the final state. 

Here, we consider the approach proposed in Ref. [5] where the double vector meson pro- 
duction has been estimated in terms of the non-forward solution of the BFKL equation [6]. It 
is based on the approach introduced in Refs. [18, 19] to describe the process 7p — > VX at 
large t, which reproduce the measured differential cross sections for longitudinal p meson and 
for J/vp production very well [20, 21]. In what follows (Section 2), we shortly present the main 
formulae of the approach proposed in Ref. [5] and its further application for the double meson 
production at UPC's. We analyze the theoretical estimates for the energies of RHIC and LHC 
and its dependence on atomic nuclei. The event rates are also computed, giving the magnitude 
of the contributions to the possible final states at the coherent collisions. Finally, in Section 3 
we discuss the main backgrounds for the photon-photon process in UPC's and summarize our 
results. 
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COLLIDER 


ION 


VS NN (GcV) 


C A a (cm^s" 1 ) 


RHIC 


Si 


250 


4.4 x 10 28 




Au 


200 


2.0 x 10 26 


LHC 


Ca 


7200 


5.2 x 10 2y 




Pb 


5500 


4.2 x 10 26 



Table 1: Luminosities and beam energies for AA collisions at RHIC and LHC. 

2 Vector meson pair production in the BFKL approach 

In the last few years several studies has analyzed the double meson production in two-photon 
processes considering different approximations for the QCD dynamics [4, 16, 7, 5, 23, 24]. In 
particular, in our previous paper in Ref. [7], we have performed a phenomenological analysis 
for the double J/^ production using the forward LLA BFKL solution. In that case, the 
hard scale was set by the charm quark mass. There, we also studied the possible effects of 
corrections at next to leading approximation (NLA) level to the BFKL kernel investigating the 
slow down of the effective hard Pomeron intercept. Afterwards, in Ref. [5] the non-forward 
solution was considered for a larger set of possible vector meson pairs, where the large t values 
provide the perturbative scale. In that paper the double vector meson production in real 
photon interactions was studied, the t-dependence of the differential cross section was analyzed 
in detail and the total cross section for different combinations of vector mesons was calculated 
using the leading order impact factors and BFKL amplitude. More recently, two other studies 
on the process 7*7* — > VV has appeared in literature [23, 24]. In the first one [23], the leading 
order BFKL amplitude for the exclusive diffractive two-p production in the forward direction is 
computed and the NLA corrections are estimated using a specific resummation of higher order 
effects. In the last paper [24], the amplitude for the forward electroproduction of two light 
vector mesons in NLA is computed. In particular, the NLA amplitude is constructed by the 
convolution of the 7* — > V impact factor and the BFKL Green's function in the MS scheme. In 
addition, a procedure to get results independent from the energy and renormalization scales has 
been investigated within NLA approximation. In what follows we shall consider the approach 
proposed in Ref. [5] , which is valid for real photon interactions and has been compared with the 
HERA data when applied for photon-proton collisions. As it is based on the leading logarithmic 
approximation (LLA) for the impact factors and BFKL amplitude, we should consider our 
results as an upper bound for the total cross sections. 

Following Ref. [5] we have that the differential cross section for the process 77 — > ViV 2 is 
given by 

6^(77 -> V X V 2 ) 167T . . , 2 , 
It = ^bfklM| , (3) 

where z = (3a s /27r) \n(W 2 / A 2 ) , r = \t\/My, My is the mass of the vector meson and A 2 is a 
characteristic scale related to M v and \t\ (See Ref. [5] for details). Moreover, in the leading 
logarithm approximation (LLA) the coupling constant is fixed and we assume a s = 0.2 in our 
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RHIC 




LHC 






Si 


Au 


Ca 


Pb 


PP 


10" 2 (0.42) 


1.2 (0.2) 


4.3 (2 x 10 3 ) 


432 (181) 


LULU 


6 x 10~ b (3 x 10~ 3 ) 


8 x 10~ 3 (10~ 3 ) 


0.03 (15) 


3 (1.2) 




2 x 10~ 3 (0.01) 


0.02 (4 x 10" 3 ) 


0.14 (74) 


14 (6) 


put 


7 x 10~ 4 (0.03) 


0.1 (0.02) 


0.3 (174) 


34 (14) 


p<p 


10~ 3 (0.06) 


0.16 (0.03) 


0.7 (378) 


73 (31) 


LU(j) 


6 x 10~ 4 (0.03) 


0.013 (3 x 10~ 3 ) 


0.06 (32) 


6(3) 



Table 2: Double light mesons: The total cross section (number of events/month) for double 
vector meson production in two-photon reactions in UPC's at RHIC and LHC energies (t min = 1 
GeV 2 ). Cross sections in units of nb. 

calculations, which is determined from a fit to the HERA data (for details see [21]). The leading 
order BFKL amplitude and lowest conformal spin (n = 0), is given by [6] 

•W(.,t) = ^J dv ^^^ e^e^Q^lf^Q^^ (4) 

where Q± is the momentum transfered, t = —Q]_, (the subscript denotes two dimensional 
transverse vectors) and 

X (v)=4XeU(l)-if>(l + uX) (5) 



,2 

is proportional to the BKFL kernel eigenvalues [22], with ip(x) being the digamma function. 

The quantities I2 Vi are given in terms of the impact factors T lV . and the BFKL eigenfunc- 
tions as follows [19], 

[Q±) ~ 1 S Q* r(l/2 + ti/)U ) A/2-,00 2vrA4M, 2 J W 



Ql r(l/2 + ii/) V 4 J Ji/2-ioo 
T 2 (l/2 + u)r(l/2 - u/2 - w//2)r(l/2 - u/2 + iv/2) 
X r(l/2 + u/2 - ii//2)r(l/2 + u/2 + ii//2) ' 

where Cj = 3r^My./o; em . Notice the non-relativistic approximation for the impact factor is 
considered here. The differential cross section can be directly calculated substituting the above 
expression in Eq. (4) and evaluating numerically the integrals. The total cross section will be 
given by 

CT ( 77 ^ Vl V 2 ) = d\t\ , (7) 

J\t\min CL\Z\ 

where \t\ m i n is the minimum momentum transfer in each particular process. In order to use 
a perturbative approach for the interaction of real photons the presence of a hard scale is 
necessary. It can be the mass of the vector meson or the momentum transfer in the process. 
Consequently, for heavy vector meson production we can assume \t\ min = 0, while for light 
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Vi Vo 


RHIC 




LHC 






Si 


Au 


Ca 


Pb 


pp 


3.2 (140) 


340 (68) 


2821 (10 b ) 


2 x 10 b (10 b ) 


IOLO 


0.02 (1) 


2 (0.5) 


20 (10 b ) 


1.6 x 10 3 (690) 


<m> 


0.02 (0.9) 


1.6 (0.3) 


20 (10 b ) 


1.6 x 10 3 (704) 


pu) 


0.5 (23) 


55 (11) 


500 (2.10 b ) 


4 x 10 b (1.6 x 10 b ) 


p(f> 


0.2 (10) 


22 (4.5) 


230 (10 b ) 


2 x 10 b (1.6 x 10 b ) 


UJ(j) 


0.02 (0.9) 


2 (0.4) 


19 (10 b ) 


1.6 x 10 3 (667) 



Table 3: Double light mesons: The total cross section (number of events/month) for double 
vector meson production in two-photon reactions in UPC's at RHIC and LHC energies (t min = 
GeV 2 ). Cross sections in units of nb. 

vector meson production we need to consider a lower limit for the momentum transfer. Here, 
we will assume \t\ min = 1 GeV 2 as in Ref. [5]. In principle, for smaller values of t non- 
perturbative contributions should dominate the total cross sections. In what follows we estimate 
the double vector meson production assuming \t\ m i n = when a heavy vector meson is present 
and \t\ m i n = 1 GeV 2 for double light vector meson production. For comparison, we also present 
our predictions for the latter process assuming \t\ min = 0, which must be used with some 
caution, since we are extrapolating the perturbative approach for very small momentum transfer 
where its direct application is questionable. However, as we expect that the t-dependence will 
not be strongly modified at small t, it can be used as a rough estimate. 

Lets now compute the double meson production in ultraperipheral heavy ion collisions 
within the BFKL approach. Using Eq. (3) into Eq. (7) and substituting it in Eq. (2) we 
can calculate the total cross section for the double vector meson production at different values 
of the center of mass energy and distinct ions. Moreover, using these results we can estimate 
the expected number of events for the RHIC and LHC luminosities considering the collision of 
different ions at distinct center of mass energies. In Table 1 we present the luminosities and 
beam energies for AA collisions at RHIC and LHC which we use in our calculations in order to 
estimate the number of events. Moreover, we assume the standard 10 6 s/month heavy ion run 
at the RHIC and LHC. Notice that the number of events by year are different for RHIC and 
LHC, since RHIC is a dedicated heavy ion mode with a run of 10 7 s/year and LHC is planned 
to have one month (10 6 s) of run in its ion mode. Moreover, we have that the center of mass 
energy at LHC is factor 20 larger than the RHIC one. 

In the Tables 2, 3, 4 and 5 we present our predictions for the production of different combi- 
nations of two vector mesons and distinct energies and ions. In particular, in the Tables 2 and 3 
we present our results for double light meson production assuming t m i n = 1 GeV 2 and t m i n = 0, 
respectively. The first aspect which should be emphasized is the large difference between the 
predictions shown in the two tables. As already shown in Ref. [5] , the differential cross sections 
peaks at small values of t, which implies a large contribution from the region of small values of 
momentum transfer. We have that independently of the choice for the minimum momentum 
transfer the production of light vector mesons is predicted to be huge, mainly at LHC energies 
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and heavy ion collisions, which become the analysis of this process feasible. In comparison 
with our previous calculations for double-p production, we have that our predictions using the 
BFKL approach are somewhat smaller or same order than obtained in Ref. [15] in terms of the 
pomeron-exchange factorization relations. However, the comparison depends on the particular 
choice for the photon parton distribution function (pdf) used in the pJ/^ calculation and on 
the treatment of the double J/ty production. Namely, in the pomeron-exchange factorization 
relations, <r(77 — > pp) oc [(7(77 — > pj /\I/)] 2 / '(7(77 — > J/^fJ/^f). The result presented here (for 
t m in = 0) is equal to prediction in Ref. [15], <r(AA — > ppAA) = 200 pb, obtained choosing the 
GRS NLO parameterization for the photon pdf and a NLO BFKL treatment for the double 
J/^f cross section. 

In Table 4 we present our predictions for the production of double light-heavy mesons. We 
have that the predicted cross sections for pJ/^ production in Pb-Pb collisions at LHC are 
in agreement with the results obtained in Ref. [15] using a distinct approach. In Ref. [15], 
the pJ '/\& cross section at photon- level was computed in the double logarithm approximation 
(DLLA) of pQCD and its magnitude depends on the particular choice for the photon pdf. The 
result presented here is similar to those calculations, a(AA — > pJ/^AA) 1 — 2 ph, when 
GRS NLO parameterization for the photon is considered. Finally, in Table 5 we present our 
predictions for the double heavy vector meson production. We have that the cross sections 
diminish for heavier vector mesons. In particular, we have that our result for the double J/\P 
production (A = Pb and LHC energy) is similar the one obtained in Ref. [7] where we have 
used a phenomenological analysis also based on the BFKL approach. We have still predicted 
a reasonable rate of production, allowing future experimental analyzes, even if the acceptance 
for the J/\I/ detection being low. On the other hand, the detection of double T production 
will be a hard task. Due to the higher integrated luminosity for light ions, the event rates are 
enhanced in this case despite the integrated cross section being small. The present estimate for 
the double J/^f production is consistent with our previous phenomenological analysis using the 
forward BFKL solution in Ref. [7], a(AA — »■ 2J/^fAA) w 60 nb. There, the main uncertainty 
to the overall normalization was the value for the diffractive slope Bj/y which is clearly 
determined in the present treatment. 

In order to evaluate how much the BFKL resummation amplifies the cross sections, we can 
compare our predictions which those obtained at Born level. We obtain that <jBom{AA — > 
ppAA) = 6 (736) nb (t min = 1 GeV 2 ), a Born (AA -> pJ/VAA) = 0.11 (129) nb (t min = GeV 2 ) 
and a Born (AA -> J/#J/#AA) = 0.001 (11) nb (t min = GeV 2 ) at RHIC (LHC) energies. It is 
important to emphasize that although the 77 — > V1V2 is energy independent at the Born level, 
there is a dependence on energy which comes from the equivalent photon flux. We can observe 
that the BFKL effects imply cross sections which are a factor 2-3 larger than the Born one, 
while at LHC this factor is > 10. 
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Vi Vo 

* 1 * z. 


RHIC 




LHC 






Si 


Au 


Ca 


Pb 




8 x l(T a (0.34) 


0.25 (0.05) 


15 (7 x 10 a ) 


1200 (504) 




6 x 1(T 4 (0.03) 


0.02 (4 x 10~ 3 ) 


1.3 (657) 


101 (42) 




4 x 10~ 3 (0.18) 


0.08 (0.01) 


2 x 10 4 (10 Y ) 


8 x 10 4 (3 x 10 4 ) 


P T 


10~ 6 (5 x 10~ 5 ) 


3 x 10~ e (6 x 10~ 7 ) 


0.02 (10) 


1.3 (0.5) 


ujT 


10~ 7 (4 x 10~ 6 ) 


2 x 10~ 7 (5 x 10~ 8 ) 


10~ a (0.85) 


0.1 (0.04) 




2 x 10~ 7 (8 x 10~ 6 ) 


4 x 10- y (10- y ) 


2 x lO" 3 (1.5) 


0.2 (0.08) 



Table 4: Double light-heavy mesons: The total cross section (number of events/month) for 
double vector meson production in two-photon reactions in UPC's at RHIC and LHC energies. 
Cross sections in units of nb. 





HEAVY ION 


J/V J/V 






RHIC 


SiSi 


2 x lO" 3 (10- 2 ) 


1.6 x lO" 7 (7 x lO" 6 ) 


8 x lO" 10 (4 x lO" 8 ) 




AuAu 


2 x 1(T 3 (5 x 1(T 4 ) 


2 x 10~ 7 (3 x 10~ 8 ) 


6 x HT 11 (10" 11 ) 


LHC 


CclCcl 


0.74 (387) 


3 x lO" 3 (1.8) 


8 x 10~ 5 (0.05) 




PbPb 


61 (26) 


0.25 (0.1) 


5 x 10~ 3 (2 x lO" 3 ) 



Table 5: Double heavy mesons: The total cross section (number of events/month) for double 
vector meson production in two-photon reactions in UPC's at RHIC and LHC energies. Cross 
sections in units of nb. 

3 Summary 

Before to present a summary of our results, lets briefly discuss the main backgrounds for the 
photon-photon processes in UPC's. An important background is given by the photonuclear 
interactions, since these reactions have similar kinematics. For the processes considered here, 
the diffractive vector meson production in photon-pomeron interactions [26, 27, 10] should 
contribute significantly. In particular, because the cross section for this process is large, as 
predicted and recently confirmed by the experimental result obtained by the STAR collabo- 
ration [28], the probability of having double (independent) production of vector mesons in a 
single nucleus-nucleus collision and associated to multiple interactions, is non-negligible [27]. 
Moreover, the two-photon cross sections have been estimated as being at least 10 3 smaller than 
the corresponding photoproduction cross sections [25], which make the experimental separa- 
tion between the two interactions very hard. Our calculations indicate that the inclusion of 
the QCD Pomeron effects implies higher cross sections at two-photon level and, consequently, 
larger cross sections in ultraperipheral collisions. Therefore, the inclusion of these effects im- 
plies that, in general, the contribution of two-photon interactions is non- negligible. Moreover, 
in principle, an analysis of the impact parameter dependence should allow to separate between 
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the two classes of reactions, since two-photon interactions can occur at a significant distance 
from both nuclei, while a photonuclear interaction must occur inside or very near a nucleus. We 
notice that the experimental separation between the two classes of processes is an important 
point, which deserves further precise studies. 

As a summary in this paper we have estimated the double vector meson production in 
ultraperipheral heavy ion collisions within the BFKL approach. We have considered the non- 
forward solution of the BFKL equation and calculated the total cross section and event rates 
for different combinations of final states, energies and heavy ions. Our results demonstrate 
that the number of events predicted at future colliders is large, allowing a future experimental 
analysis of this process. It is important to emphasize that in our calculations we have considered 
the leading order solution of the BFKL equation as well as leading order impact factors. In 
principle, the inclusion of the next-to-leading order corrections implies a smaller growth of the 
cross section with the energy. Similarly, saturation effects can contribute significantly, reducing 
the total cross section. Consequently, our predictions should be considered as an upper bound. 
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